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Abstract —In this paper we develop a new technique for 
estimating fine clock errors and range between two nodes 
simultaneously by two-way time-of-arrival measurements us¬ 
ing impulse-radio ultra-wideband signals. Estimators for clock 
parameters and the range are proposed that are robust with 
respect to outliers. They are analyzed numerically and by means 
of experimental measurement campaigns. The technique and 
derived estimators achieve accuracies below 1 Hz for frequency 
estimation, below 1 ns for phase estimation and 20 cm for 
range estimation, at 4 m distance using 100 MHz clocks at both 
nodes. Therefore, we show that the proposed joint approach is 
practical and can simultaneously provide clock synchronization 
and positioning in an experimental system. 

Keywords: clock synchronization, range estimation, robust esti¬ 
mators, ultra-wideband 

I. Introduction 

A common time reference among nodes along with node 
positioning enables coordination in time and space domains. 
By such coordination the nodes can execute tasks efficiently in 
internet-of-things (IOT), machine-to-machine (M2M) commu¬ 
nication, or similar scenarios 0-0- In certain scenarios time 
synchronization among nodes is necessary for positioning. 
Specifically, it is needed in positioning techniques where time 
of arrival (TOA) measurements are required and measurements 
are performed using a clock. As an example, in the commer¬ 
cially available Ubisense system, sensor nodes are synchro¬ 
nized using cable in order to perform time difference of arrival 
(TDOA) measurements of ultra-wideband (UWB) pulses from 
tag nodes 0. In these types of systems, replacing the cabling 
with wireless synchronization will reduce the installation cost 
in stationary environments and increase in flexibility will open 
up new mobile applications. Positioning and effects of clock 
parameter mismatch have been addressed jointly before in (0- 
0 and the references therein. In 0 , optimal algorithms were 
derived for joint positioning and synchronization with anchor 
uncertainty in position and time. In [6|, the authors proposed 
distributed algorithms for positioning and synchronization in 
UWB ad-hoc networks. In (7], the effect of clock skew was 
discussed and estimators for TDOA-based positioning were 
formulated. A few recent works on joint ranging and clock 
synchronization are mentioned in (8|-(TT). In Q, authors have 
proposed a network wide clock synchronization and ranging 
mechanism by pairwise timestamp measurements aided with 
passive listening. In fl0| , authors have proposed clock syn¬ 
chronization and localization along with a low cost testbed 
for its demonstration. 

Our main contribution in this paper lies in suggesting a mea¬ 
surement mechanism which provides information about clock 
parameters and range between nodes. Further, we propose joint 


estimators to precisely estimate range and clock parameters 
from the measurements. The estimated range can subsequently 
be used for positioning using standard techniques (T0-(T4). 
Such an approach would simplify physical infrastructure and 
obviate the need for e.g. laying cables for synchronization in 
positioning systems. 

Similarly, wireless clock synchronization can be useful in 
several other applications; which either presently rely on wired 
synchronization or work asynchronously but exhibit perfor¬ 
mance loss in the absence of synchronization. For instance, 
many packet-based radio systems operate asynchronously and 
suffer packet collisions. Joint positioning and time synchro¬ 
nization can help reduce such losses in a network. Tradi¬ 
tionally, time synchronization is a necessity for a variety 
of networks. Network time protocol (NTP) is a well-known 
mechanism to synchronize clocks over the internet G3- 
Similarly, time synchronization protocol for sensor networks 
(TPSN) is in use (16). In the reference broadcast synchro¬ 
nization (RBS) protocol, time synchronization is achieved by 
broadcast beacon and receiver timestamp exchanges p7) . In 
all these network synchronization mechanisms, the clocks are 
eventually synchronized between two nodes in a network. In 
this paper, by contrast, we avoid timestamp exchanges, and 
hence communication overhead, using round-trip time (RTT) 
measurements. This enables simultaneous estimation of the 
clock parameters and range between two nodes. 

To achieve high accuracy in wireless clock synchronization 
it is necessary to estimate the delay arising from the signal 
propagation time between nodes. Here we consider using 
UWB as the wireless technology and demonstrate results using 
a UWB-based radio with a very accurate time measuring 
device or chronometer. Our previous work has demonstrated 
time measurement precision up to 50 ps using a time-to- 
digital converter (TDC) |18;]. Furthermore, usage of UWB in 
very accurate range and position estimation is already well 
known §, jig, 00, 03 . Under line of sight conditions 
the signal propagation delay is directly proportional to the 
distance between the nodes. Estimating this delay therefore 
enables ranging as a by product of clock synchronization. 

Very few clock parameter estimation results have been 
verified experimentally in the literature. In [|T0, measurements 
were captured using Berkeley motes and logic analyzers and 
verified offline by synchronization algorithms. The timescale 
and accuracy in these experiments are in microseconds. 
Whereas, in our proposed system the timescale and accuracy is 
in the order of a nano-second, as shown in subsequent sections. 
In (20) , a timestamping system for temporal information 
dissemination is proposed achieving sub-nanosecond accuracy 


with UWB pulses. The timestamping system is integrated in 
a ZigBee platform for wireless sensor networks. However, in 
[20), one-way time-of-arrival (TOA) is considered, and range 
estimation is not supported, differently from our proposed 
technique. 

In the literature various theoretical models and methods 
have been discussed for clock synchronization over wire¬ 
less networks, cf. 121|-|23|. In [1221, authors have suggested 


message-passing methods for network clock estimation. The 
relative clock model between two clocks adopted in next 

the 


section is discussed in |23| in greater detail. In 
Cramer-Rao bounds were derived on the parameters for this 
clock model. In subsequent sections we will propose the usage 
of a clocked delay in the slave node. In our previous work, 
the generated delay is assumed to be fixed and analog G3- 
[ [25) , (26) . Here, by contrast, we turn the clock into a delay¬ 
generating device producing informative measurements. 

Since the proposed technique estimates the range between 
a master and a slave node, it can be used as a fundamental 
building block of scalable and cooperative distance-based 
techniques for network positioning. Specifically, an example 
of a widely-used technique exploiting node-to-node ranging is 
based on multi-dimensional scaling (MDS) (27). Furthermore, 
cooperative message-passing methods based on Bayesian net¬ 
works have been developed in (12) , and experimentally eval¬ 
uated on a UWB ranging platform in (28) . An analysis and 
implementation of such methods is beyond of the scope of the 
present paper. 

The remainder of the paper is organized as follows: Section 
|H| describes the considered system model for the paper. In 
section II-B the clock parameter measurement mechanism 


is introduced. This section also introduces the measurement 


model. In section III estimators are proposed for estimating 
clock parameters and range from the measurements. Section 
[Tv] is composed of a numerical evaluation of the model 
proposed in section [II] Section [V] provides results from exper¬ 
imental setup and performance of estimators on the acquired 
data. In section VI we conclude the paper. 

Notation: 1 denotes a column vector of Is. a © b is 
the element-wise, or Hadamard, product between vectors a 
and b. x 1 / 2 is the element-wise square-root of vector x. 
||x||w = Vx T Wx where W is a positive semidefinite matrix. 
The modulus function is written compactly as mod a (x) = 
modulo(x, y). 


II. System Model and Problem Formulation 

In this section, we first describe the principle of operation 
of the considered system, then we illustrate the proposed clock 
measurement mechanism. Subsequently, we define the related 
signal model and present a basic estimator of the parameters 
of interest, based on existing literature. 



Fig. 1. Principle of operation of the considered system: the master node sends 
a PING signal and the slave node responds. The master node is equipped with 
a time-to-digital converter capable of measuring the round-trip-time. 


PlNG signal and the reception of a RESPOND signal, sent by 
the slave, after a predetermined delay. A conceptual overview 
is given in Fig. [1] This time interval is denoted as the round- 
trip-time (RTT) which we express as a function y of the 
unknown range between the master and slave p and the delay 

15 - , 2 

y = S+-p, (1) 

c 

where c is the speed of light in meters per second. Convention¬ 
ally, the delay is an ideal analog implementation and known 
at the master, cf. on, m- Here, however, we assume that 
all timed events are based on local clocks at each node. As 
we will show below, this enables the joint estimation of both 
range and parameters for clock synchronization. 

Time is measured using a clock by counting a number of 
periods of the clock signal 

„ n 
Co — — + p 

Jo 

where n is the number of clock cycles, fo is the nominal 
clock frequency and P G [0,1//o) is the initial phase of the 
clock when the measurement begins. Any fractional deviation 
from the nominal clock frequency is termed as skew of the 
clock, represented usually as a > 0. Time is measured using 
a skewed clock, expressed in the nominal frequency fo as 

77 

C = a y o +V’ 

and the actual frequency / deviates from nominal clock 
frequency given by a = fo/f- Equivalently the above model 
is widely used as Q, (29), (30): 

C = at + <p, 

For simplicity and to establish the notation, let us first 
consider a standard linear clock model of the observed times 
at master and slave nodes. 


A. Principle of Operation: Round-Trip Time 

We consider two nodes equipped with a wireless transmitter 
and a receiver. One being called the master and the other 
called slave | 181. We assume that the master node is capable 
of measuring the time interval between the transmission of a 


and 
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C s = a s t + p s , 

where a m and a s denote their clock skew with respect to 








































Fig. 2. Relative clock error models, (C m ,C s ) in J 2 J a s /a.m denotes the 
clock skew and O initial phase in radians. 


nominal frequency f 0 and ip is the initial phase offset. 


The clock parameters are assumed time-invariant 124]. If 
f m and f s denote the clock frequencies at the master and 
slave, then we define the frequency difference fd = fm — fs- 
The relative clock skews will then be parameterized by fd as 
Us fam = fm/fs = fm/{fm - fd)- The time observed at the 
slave C s can then be expressed relative to a reference, i.e. the 
master node clock C m , 
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where the second term is the relative clock offset in seconds 
and we correspondingly define the relative clock offset in 
radians <j> as 

4> = 2n(f m - fd) (<Ps- — e [0, 2tt). 

\ Q-m J 

This model is illustrated in Fig. [2] Given that C m is the 
reference clock, f m is known and to perform synchronization 
of master and slave clocks, only estimation of fd and cf> is 


necessary 1311. 

If we assume that the slave generates the delay by counting 
a fixed and known number of clock cycles, equivalent to 
quantization in the time domain, the nominal delay So becomes 
dependent on the frequency difference and can be written as 

l{C' s -C s ){f m -f d )\ 

O0 (Jd) — -7 ——? - j 

Jm Jd 

where C s and C' s denote the time of reception of the PlNG 
and transmission of the RESPOND signal at the slave. 

On the other hand, the master is equipped with a TDC, 
and is therefore able to measure time intervals with a much 
higher resolution than its clock period. Based on this mecha- 


Fig. 3. A space-time diagram depicting the clock edges (vertical arrows) at 
the slave (top row) and master (bottom row). The diagram relates the quantities 
in {TJ and j3j. Specifically, the remainder term //.(/.) in j3j is highlighted. 


nism, which we will describe in more detail below, the time 
measured at the master is given by an integer number of slave 
clock cycles (i.e. <5o) plus a remainder term, which is naturally 
modeled by a periodic modulus operation. In the following 
subsection we present a model of this term, denoted h(t ; fd, <j>) 
which is dependent on the unknown clock parameters fd and 
<j>. Therefore, we can write the delay as 

S = S 0 {fd) + h{t-f d ,(j)) (3) 

and in combination with ([T]( the round-trip-time measurements 
will be informative of the range p to the slave but also of the 
relative clock parameters of the slave, fd and <f>, via 6. 


B. Clock Measurement Mechanism 


Assuming the frequency difference is small, i.e. fd <C f m , 
we can neglect the dependence of <5o on fd since /, = /,„ 
fd — fm an d thus we can approximate So as a known constant, 

_ [(^ -Cs)f,„] 
do- -7-■ 

J m 

Similarly, ([3]) can be approximated at a specific time instant t 
as 


S = 


L(c' - c s )(f m - f d )j 
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mod 2n{fdt + <t>) 
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where periodic modulus operation represents the remainder 
and is dependent on the clock offset and frequency difference. 
We now illustrate how the remainder term appears at the 
master node. Consider the following example based on the 
events described in Fig. |T| 

1) Master sends a PlNG pulse on a positive edge of clock. 




































TDC receives the Start signal instantly and starts 
measuring time. 

2) On receiving the pulse, slave starts the delay generation 
from subsequent positive edge of its own clock. 

3) Slave sends RESPONSE after the delay (Jo (say, two clock 
cycles in this example) 

4) Master receives the Response, TDC gets the Stop 
signal and measures the RTT. 

This sequence of events is also illustrated in the space-time 
diagram of Fig. [5] where the remainder term h(t) is high¬ 
lighted. The master repeats the above transmission procedure 
with a fixed periodicity. The result is illustrated in Fig. [4] which 
shows the role of the clock in repeated pings. At instance 
‘5’ Ping from master is received by slave at a clock edge 
and hence results in smallest RTT. Whereas, at instance ‘6’ 
slave receives PlNG narrowly after the clock edge and hence 
waits nearly a full clock duration to start its delay count of a 
clock cycle. This results in largest RTT measurement at master. 
As can be seen, the clock offset and frequency results in a 
perodic time-varying remainder term h(t) with a sawtooth-like 
waveform corresponding to the modulus function that follows 
from rounding. 

In addition to the effect of the rounding operation, the 
remainder term in 0 is also subject to zero-mean random 
clock jitter, which we denote v. At measurement time instant 
ti we therefore model the remainder as 

h{ti\ f d , 4>) = ^mod 2vr (/ d f + cj> + v{ti)), (5) 

Z7T 

where T m = 1/ f m . Finally, inserting 0 into 0 under f d < 
fm, Q results in the following the RTT measurement model 
at time tp. 

T 2 

y{ti ) = ^mod 27r (/ d U +$ + v(ti )) + S 0 + -p + n(ti), (6) 

Z7T C 

where n(£j) denotes zero-mean measurement noise. Figure [0 
shows a measurement capture of above phenomena on our 
flexible UWB radio test-bed m ■ In Section [V-E| we validate 
the model 0 by means of residual analysis, justifying the 
approximations used above for the considered experimental 
setup. 

C. Signal Model and Measurements 

Let y = [y{ti) ■ ■ ■ y(t N )\ T £ R N and t = [fi • • • t N ] T £ 
R ;V be the vector of N measurements and time instances, 
respectively. Then model 0 can be written compactly in 
vector form: 

y = h(/ d ,0,v) + <J 0 1 + — 1 + n, (7) 

C 

where 

Hfd, <P,v) = ^mod 27r (/ d t + (f)l + v). 

Z7T 

The vectors v and n contain the random noise terms with 
unknown variances and er^, respectively. The goal is to 
estimate p, f d and </> from y. 

Central to RTT measurement procedures is the device called 
time-to-digital converter (TDC) as is shown in Fig. [T] Com¬ 
mercially available TDCs can measure time with resolution 



sample index [i] 


Fig. 5. Experimental data record showing the “sawtooth” behavior of the 
RTT measurements. The data was obtained with the master and slave antennas 
placed 0.5 m apart, <5 q = 4900 ns. 


up to 50ps |32j. As shown in Fig. [I] TDC measures the time 
difference between two events and thus measures RTT. The 


TDC has its own clock for coarse counting 1321. 

Here we wish to highlight a few important features: 


• Measurement model 0 arises from sampling the slave 
clock over a wireless medium and measuring the relative 
clock drift using a TDC at master node. 

• No timestamping and hence no data transfer is required 
between master and slave devices. 

• Traditionally, clock parameters are estimated as shown 
in Fig. [2] where measured time monotonically increases. 
Long time measurements using a clock require high 
accuracy. Whereas, RTT measurements as seen in Fig. [5] 
require short time measurements, hence clock accuracy 
and stability is a lesser concern. When measuring smaller 
time duration the TDC clock skew is assumed to be 
negligible. 

• Clock parameter information is available in the time 
domain as well as the frequency domain. The latter fact 
facilitates extracting information using frequency-based 
methods. 


Next, we introduce a conventional estimator for estimating 
fd, 0 and p by unwrapping the saw-tooth measurement 0; 
an early work is J29|, cf. discussion in ]24), (33). 


D. Unwrapped Least Squares (ULS) 

We consider the model in 0 and formulate a three-step 
unwrapped least-squares method (ULS). First, by approximat¬ 
ing the arithmetic average of the modulus term in 0 as the 
amplitude of the sawtooth-like waveform, T m / 2, we obtain 

2—1 
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Fig. 4. A diagram depicting the clock mechanism. Note that the RTT measurements, plotted below, vary depending on the shift of the slave clock edge 
relative to the master clock. 


From this we form an approximate least-square estimate of 
the range 
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p = arg mm 
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_ T m 2 p 


L j- 

= o \y ~ -w - 


Then, we normalize the measurements y(ti) to the [—7 t,7t)- 
interval as follows: 

Z(U) 4 (y( U ) - y) . 

-f m 

Next, we apply a ‘phase-unwrapping’ algorithm, denoted by 
U{ z) where z = [z{h) ■ ■ ■ z{t N )] r §4j, (29], |34|. This 
algorithm adds multiples of ±27 t rad to the phase when the ab¬ 
solute difference between consecutive elements is greater than 
7 r rad. Finally, we obtain & and f d by solving the following 
linear least-squares problem on the unwrapped samples: 

[fd,<t>] = argmin \\U(z) - f d t - <j)l\\ 2 2 . (g) 


Whilst the ULS estimator is computationally efficient, erro¬ 
neously unwrapped samples by U{-) which occur in poor 
signal conditions lead to performance degradation. 

Furthermore, as formulated it is not robust with respect to 
noise outliers and its estimates may deteriorate significantly. 
The presence of outliers is a common issue in practical 
wireless measurement systems. As an evaluation example, the 
wireless localization scenario is considered in (35). In such 
a scenario, outliers in the time-of-arrival measurements are 
present due to propagation channel conditions, and robust re¬ 
gression techniques are employed for the purpose of reducing 
the adverse impact of outliers on localization performance. 


III. Proposed Estimators 
In view of the limitations of the unwrapped least squares 
approach described above, we develop two alternative methods 
to estimate the offset frequency f d , offset phase </> and range 

p in Q. 


A. Periodogram and Correlation Peaks (PCP) 

In the first approach we exploit the periodicity of the 
signal model (6) to obtain a frequency estimate. Using the 


periodogram [36j] we can compute an estimate up to an 
unknown sign. 


fd = arg max 

fd&F 


N 

Y,y(U)e~ Mu 


where F C [0, / max ] is a uniform grid of frequencies. Under 
uniform sampling, the estimate can be computed efficiently 
using the fast Fourier transform. Note that f d determines the 
slope of the periodic sawtooth waveform in (6), cf. Fig. [5] 
The sign of f d can be resolved jointly with phase estimation 
by correlating the observed signal with sawtooth waveforms 
based on f d and —f d , corresponding to positive and negative 
slopes. 

Then using the estimate f d a waveform p(<j>, s) = 
mod^^f-s./j/t + (j> 1) is generated, for a fixed frequency sign 
s = {—1,1} and nominal phase <fi. The phase estimate, and 
the resolved sign of the frequency, is then obtained by the 
correlation peak, 

[<M] = argmax |y T p(<j>, s)|, 
se{-i,i} 


where <I> C [0, 27 t) is a uniform grid of phases. The frequency 
estimate is f d = sf d . Finally, the range estimate is obtained 
by a direct least-squares fit using 0 which results in a 
particularly simple form. Using vector notation, it can be 
written as 


p = arg min 

p 


y- 


-^POM) - <5ol- 



s) - So 1 
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We summarize this low-complexity three-step approach as 
computing the periodogram and correlation peaks (PCP). This 
approach circumvents the need for phase unwrapping the 
data and utilizes frequency domain characteristics of measure¬ 
ments. 


B. Weighted Least Squares (WLS) 

In the second approach we address the issue of robust 
estimation in the presence of noise outliers, using a weighted 

































squared-error criterion 

J{fd,<i>,p) = ||y - Hfd,<t>, 0 ) - 5 0 1 - 2c"Vl||w > ( 9 ) 


and find the minimizing arguments. Here W = diag(w) y 
0 £ R^ xN is a diagonal weight matrix. The weights w can 
be chosen to mitigate the effect of outliers by downweighting 
the samples likely to be outliers, as discussed in the next 
subsection. This makes WLS robust with respect to outliers. 
Setting uniform weights w oc 1 reduces 0 to a standard 
least-squares criterion. 

We begin by minimizing ([9]) with respect to p. This gives 

„ cw T (y - h(/ rf ,/>,0) - 5 0 1) 


Inserting the estimate © back into 0 and defining 
r(/d, <t>) = y - h(/d, <j>, 0) - (5 0 1, we can write the mini¬ 
mization with respect to fd and 0 as a two-dimensional grid 
search 


[fd,<t>] = argmin ||w 1/2 © r(/ d , 0)||1 

fd£F, $ 


|w T r(/d, (f>)\ 2 

l T w 

( 11 ) 


where F C [—/max, /max] and $ C [0, 2n) denote the grids. 
See the appendix for the derivation. 


C. Choosing the Weights for Outlier-prone Data 


For applying the WLS method derived in Section III-B 


to 

outlier prone data, we propose to choose the /-th element w t 
of the weight vector w as follows: 


Wi 


1 if \y(U) - pi /2 (y)| < 3<j mad (y) 

0 otherwise. 


( 12 ) 


where fi/ 2 (•) denotes the sample median operator, and <7 mac j 
denotes the normalized Median Absolute Deviation (nMAD), 
which is defined, cf. [35), as 


o-mad (y) = 1-483 • p 1/2 (y - lp 1/2 (y)) . 

As the signal model 0 has no trend and a constant offset, 
this choice of weights ([12) takes large deviations from the 
median to be an indicator of a likely outlier and downweights 
the corresponding samples. 

We note that other outlier-rejection strategies are possible 
but are not presented here; see [38j| for an extensive 
survey. An advantage of the proposed strategy is that it is 
automatic, since it does not require any input or parameter 
tuning by the user. Moreover, in the considered context, 
another desirable feature of the strategy is its tight integration 
in the WLS framework. Finally, the proposed strategy gives 
good results in a practical outlier-prone scenario, as we show 
in the experimental evaluation of Section [V] 


IV. Numerical Results 

In this section, we explore the properties of the model in 
([7) and evaluate the performance of the proposed estimators 
in realistic scenarios. 


A. Simulation Setup 

The numerical simulations have been performed, unless 
otherwise indicated, using the following fixed values of the 
parameters of interest: fd = —32 Hz, p = 2 m. For each 
Monte Carlo iteration, the initial phase f is set to a random 
value according to a uniform distribution in the [0, 27t) interval. 
Unless otherwise noted, a record comprised of N = 100 
samples is generated for each Monte Carlo iteration according 
to the model in 0. The true value of the master clock 
frequency is T m = 10 -8 s, the value of <5o is 5 x 10 -6 s and 
a measurement update period of T s = 10 -3 s has been used. 
Taking into account the model in (|7), we define 

SNR C ^ 10 log 10 % , SNR, ^ 10 log 10 . 

®n 

SNR C is signal-to-noise ratio over the wireless channel and 
SNR, is signal to noise ratio of the clock jitter noise |39) . 
The channel noise for simulations is generated in seconds and 
the jitter noise is the dimensionless phase noise as can be seen 
from 0. The performances are evaluated using the root mean- 
square error (RMSE) and computed by averaging over 1000 
Monte Carlo iterations. 


B. Simulation Results 

A complete characterization of the Monte Carlo simulation 
results is provided in Fig. [6] Due to the fixed resolution of 
the search grid, WLS and PCP exhibit quantization effects, 
as shown e.g. in the asymptotic constant behavior at high 
SNR in Fig. [6]a. Potentially, these effects can be overcome by 
refining estimates using standard interpolation, zero padding, 
or adaptive grid search methods. However, their performance 
is still satisfactory. In fact, all considered estimators are able 
to obtain accuracy of the order of 1 Hz or better, when SNR C 
and SNR, are greater than 20 dB. 

Furthermore, we studied the behavior of the proposed meth¬ 
ods as the number of samples N varies. The simulation results 
are shown in Fig. [6]g -[6]i. The proposed WLS estimator shows 
improved accuracy in small-sample conditions with respect 
to ULS and PCP, which is advantageous in scenarios with 
slowly time-varying parameters. The performance of ULS at 
high SNR and for large N is generally better than PCP and 
WLS, as shown in Fig. [6]a, d, and g. This is consistent with the 
expected behavior. As noise vanishes, in fact, the unwrapping 
procedure results in a linear trend, which can be estimated 
statistically efficiently by the linear least-squares method. 

Finally, in order to analyze the effect of outliers, we 
performed a set of Monte Carlo simulations substituting a 
fraction of the data with outliers, uniformly distributed in the 
[3.5 • 10 -6 ,4.9 ■ 10 -6 ] s interval. The outliers are inserted at 
randomly-selected locations, yielding the results in Fig. [6jj - 
0- The proposed WLS estimator has considerably improved 
robustness with respect to the other considered methods. More¬ 
over, it exhibits a threshold behavior, reaching a breakdown 
point when the fraction of outliers is approximately equal to 
40%. On the other hand, ULS shows a degraded performance 
even at a low fraction of outliers, in Fig. [6]g and [6]h, therefore 
validating the proposed WLS approach. PCP operates first in 
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Fig. 6. Simulation results, averaged over 1000 Monte Carlo iterations, (a)-(c) varying SNR C , with SNRj = 40dB. (d)-(f) varying SNRj, with SNR C = 30dB. 
(g)-(i) varying the number of samples N, with SNRy = SNR C = 40 dB. (j)-(l) in the presence of outliers, varying the fraction of randomly-inserted outliers, 
with SNRy = SNR C = 40 dB. 








































































































































































































































































































































the frequency domain which gives it a degree of insensitivity 
to outliers. 

C. Figure of Merit 

The numerical results in Fig. [6] show that, using the consid¬ 
ered RTT clock measurement mechanism and the proposed 
WLS estimator, 100 samples are sufficient to achieve an 
estimation accuracy of Hertz-order for clock frequency dif¬ 
ference, nano-second order for initial phase difference and 
decimeter-order for ranging in practical noisy and outlier- 
prone conditions. This can be verified from the curves in 
Fig. [6](g)-(i), where an SNR C =SNR ; = 40 dB is considered. 
Then, from Fig. |6](a)-(f) it is possible to observe that such 
accuracy level is still achieved for SNR C as low as 10 dB and 
SNRj as low as 20dB. More importantly, this accuracy is 
achieved by the proposed WLS estimator when up to 30% 
of the measured data are outliers, as shown in Fig. |6](j)-(l). 
Notice that N = 100 samples correspond to an observation 
time of 0.1 s for the considered update period of T s = 10“ :1 s, 
which we have chosen to be of the same order of magnitude 
as our practical experimental setup. Such an observation time 
is feasible for most practical applications of joint wireless 
synchronization and ranging. 

D. Sensitivity to Parameter Variations 

To evaluate the performance of the considered estimators 
when the parameter fd varies, we provide the simulation 
results in Fig. [7] showing that WLS is insensitive to the 
value of fd. On the other hand, ULS and PCP show con¬ 
siderable degradation of the frequency difference estimation 
performance as fd increases, see Fig. |7Ja). As far as the range 
is concerned, we note that the model in 0 is linear in the 
range p, thus it is expected from the theory that the behavior 
does not change significantly for different values of p. This 
is confirmed numerically in Fig. |7Jc), where the observed 
performance is approximately constant. 

V. Experimental Results 

We applied the proposed estimators to experimental data 
obtained from an in-house test-bed which is based on the IR- 
UWB technology gg. In the following, we provide a descrip¬ 
tion of the experimental setup, we evaluate the performance 
of the considered estimators (ULS, PCP, WLS) on measured 
data, and eventually we discuss the validation of the proposed 
model in 0 - 

A. Experimental Setup 

Experimental setup is shown in Fig. [8] Each of the two 
nodes was running on its own clock which was provided by a 
local oscillator in the ML505 evaluation board using a Xilinx 
Virtex-5 FPGA (40) . Since the two local oscillators were not 
synchronized, a non-zero frequency difference between the 
two oscillators was present. In the experimental setup, the 
master node was programmed to repeat the RTT measurement 
procedure with a fixed sampling rate of 5 kHz, and then 
transfer each result from the TDC to the host computer, where 


it was stored and used for offline processing. The apparatus 
and components used in the experiments are relatively low 
cost and easily available. This brings the proposed idea closer 
to practice. 

As can be seen from Fig. [8] while capturing measurements, 
we capture ground-truth values of the parameters. The ground- 
truth values are used for analyzing RMSE of the estimates 
over measured data. The ground truth of range (p) is mea¬ 
sured using a handheld laser distance meter ED- Ground- 
truth values of master and slave clock frequencies ( f m , fs ) 
are captured using Agilent 53230A frequency counter (42) . 
Ground-truth of phase (0) is obtained by capturing the clock 
state using primitive delays in FPGA. The clock of slave node 
is propagated through these delays and the received first PlNG 
in a sequence PlNGs from master triggers capture of the clock 
phase (43). 

We placed the master and slave nodes on two carts in 
line-of-sight conditions in the office-like indoor environment 
depicted in Fig. [9] For each considered distance a record of 
65356 RTT measurement values was acquired. In the experi¬ 
mental conditions, we observed an SNR C varying from approx¬ 
imately 14 dB at a distance of 0.5 m to approximately OdB at 
the maximum operating distance of about 4 m. Furthtermore, 
the SNR , is approximately 40 dB, and the fraction of outliers 
varies from about 5% up to approximately 20%. 

The amplitude of the observed sawtooth waveform is ap¬ 
proximately 10 ns as can be seen in Fig. [5] which is consistent 
with the nominal clock frequency of 100 MHz of the FPGA 
boards used. Furthermore, the period is approximately 30 
samples which, considering the sampling frequency of 5 kHz, 
corresponds to a frequency difference in the order of 200 Hz. 
Such a difference is consistent with the tolerance specification 
of the local oscillators used in the experiments. 

The raw RTT data from the experimental measurement 
platform was directly used as the input for the WLS estimation 
algorithm. Whereas, the raw data was pre-processed by a 
coarse outlier-rejection method before being input into the 
ULS and PCP algorithms. In particular, the pre-processing 
outlier rejection method consists of the identification of the 
outliers by means of the same criterion in ( fl2| ), and in the 
following substitution rule: if the outlier is isolated, namely the 
preceding and following samples are not outliers, then a linear 
interpolation is performed. Specifically, the outlier sample is 
replaced by the average of the following and preceding sample. 
Conversely, if the outlier is not isolated, i.e. if the preceding 
or following samples are also outliers, the outlier sample is 
replaced by the median of the entire record. 

For validation purposes, we also performed a characteriza¬ 
tion of ULS and PCP without using the pre-processing step. 
Results of such characterization, not presented here for brevity, 
show a considerable performance degradation of one order of 
magnitude or more with respect to WLS. 

Therefore, in the following subsections, we provide results 
from experimental tests where the pre-processing method is 
applied to ULS and PCP. 



Fig. 7. Numerical simulations: RMSE performance of the proposed estimators for different values of the parameter /,/. Here, SNR C = 20 dB. SNRj = 20 dB. 
N = 200, and no outliers are present. Results are obtained by averaging over 1000 Monte Carlo iterations. 


B. Frequency Difference Experimental Results 

We acquired RTT data from the system at four master- 
slave distances from lm to 4 m at lm steps. Simultaneously, 
the related ground-truth data, i.e. the true master-slave clock 
frequency difference, was acquired using the setup above 
and found to be approximately —30 Hz with small variations 
among datasets at different distances. 

Subsequently, the acquired data was segmented into records, 
each consisting of 1000 consecutive samples. Each record was 
processed using the three considered methods. The resulting 
estimates were compared with the ground-truth, to obtain the 
results shown in Fig. El The proposed WLS method outper¬ 
forms both the other considered methods, and its performance 
is relatively insensitive to variations in range. Moreover, the 
global RMSE, computed taking into account data acquired 
at all distances, is equal to 0.96 Hz for the proposed WLS 
estimator. Thus, the capability of sub-Hertz relative frequency 
measurement in experimental conditions is demonstrated. 

Furthermore, a data record illustrating the presence of 
outliers is shown in Fig. 10 The presence of outliers can be 
explained by analyzing the behavior of the receiver employed. 
In our experimental setup, in fact, the energy-detection receiver 
cannot distinguish between a proper response signal and an 
external RF interference signal. Therefore, when external in¬ 
terference bursts cross the pre-defined threshold in the energy 
detector, they cause spurious detections. Furthermore, the 
wideband noise which is inherently present in UWB systems 
can also cause spurious detections. Such spurious detections 
cause outliers during the RTT measurement procedure. 



Fig. 8. Diagram of the experimental setup. 



Fig. 9. Picture of the experimental setup. 


C. Phase Experimental Results 

Figure [12] shows error between measured phase and esti¬ 
mated phase for various estimators. Since the phase parameter 
is same for one set of measurements, computing RMSE of 
phase estimate would require many set of measurements which 
is practically difficult and time consuming. Hence, instead of 
RMSE plots, a few phase error estimates are shown in the 
figure. As can be seen, the phase estimation error is of the 
order of 1 ns for WLS and PCP estimators. As expected from 
simulation results, phase estimation error of ULS is very large. 


D. Range Experimental Results 

In order to relate the measured RTT to the master-slave 
range, we performed a calibration procedure based on that 
described in detail in G3- The calibration procedure removes 
various systematic delays, particularly delay in RF front end 
of the transceivers. In particular, 1000 RTT samples were 
acquired at a set of known distances from 0.5 m to 4.5 m 
at 0.5 m steps. Then, the calibration curve was calculated by 
applying a fifth-order polynomial fitting to the experimental 
data. The range estimate characterization was then performed 













































































Fig. 12. Phase error measurements. 


Fig. 10. Experimental data record showing outliers denoted with magenta 
circular markers. The data was obtained with the master and slave antennas 
placed 5 m apart. The dashed green line denotes the sample mean, which is 
clearly influenced by the outliers, and the dotted red line denotes the sample 
median which is robust against outliers. 
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Fig. 11. Experimental results for the estimation of frequency difference, (a) 
RMSE of the frequency difference estimates vs range, for the three considered 
methods, (b) box plot of the frequency deviation estimation error vs range, 
for the proposed WLS estimator. The blue edges of each box denote the 25th 
and 75th percentiles, the central red mark denotes the median, and the top and 
bottom ends of the black dashed lines denote the maximum and minimum 
values, respectively. 
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Fig. 13. Experimental results for the estimation of range, (a) RMSE of the 
range estimation error vs actual range, for the three considered methods, (b) 
box plot of the range estimation error vs actual range, for the proposed WLS 
estimator. The blue edges of each box denote the 25th and 75th percentiles, 
the central red mark denotes the median, and the top and bottom ends of the 
black dashed lines denote the maximum and minimum values, respectively. 


on an independent dataset, acquired at four distances from 
lm to 4m at lm steps. We denote the latter dataset as the 
validation dataset. For each distance in the validation dataset, 
we acquired multiple 1000-sample records and processed each 
record using all the considered methods. The range estimation 
results are summarized in Fig. [13] As clearly noticeable 
from Fig. |~i~3] i, the difference between the three considered 
estimation methods is less than 1 cm and, therefore, negligible. 
The resulting global RMSE, which is calculated over the entire 
validation dataset, is 17 cm for the proposed WLS method. 
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Fig. 14. Autocorrelation function over the residuals of WLS, with 99% 
confidence interval of a white noise sequence. Measurements at 4 m. 


E. Residual Analysis 

Figure [14] shows the estimated autocorrelation function over 
the residuals of the WLS estimator. The residuals are obtained 
after subtracting noiseless signal template as in 0 using WLS 
estimates from the measured data at 4 m. The procedure of 
analyzing residuals is similar to as discussed in |44) . Fig. 14 


also plots the 99% confidence bounds for a white noise 
sequence. As can be seen, virtually all ACF points fall within 
these bounds except at the zero lag which thus validates the 
model proposed in 0- 


VI. Conclusion 

In this paper we have proposed a new method to generate 
wireless measurements for joint clock parameter estimation 
and range estimation between the clocks of two nodes. We 
have formulated a mathematical model which captures the 
measurement generating phenomena. Further, we have devel¬ 
oped two robust estimators (PCP and WLS) with different 
complexities and have compared their performance with the 
popular ULS estimator. The estimators were compared by 
simulations using the proposed model. We suggested outlier 
rejection techniques to handle real-world measurements and 
compared the estimator performances using data generated by 
a UWB testbed, reporting favorable experimental RMSE. 

We envision that this method has a two-fold use: firstly as a 
building-block in positioning systems which require wireless 
node-to-node synchronization. Secondly, this method will give 
rise to applications where accurate wireless clock synchroniza¬ 
tion would be beneficial, such as in the internet-of-things, and 
the machine-to-machine communication frameworks. 


Appendix 

To derive the WLS estimator in the efficient form of ©■ 
we first define the matrix 

n = i(i T wi)" 1 i T w = 

1 w 


where the equality follows from W1 = w. Next, recall that 
r(/d, <j>) = y - h(f d , </>, 0) - <5 0 1. Inserting ([T0]> back into ]9]>, 
yields a concentrated cost function which we can re-write: 

J(fdA,p) = ll r -nr||^ 

= r T Wr - 2r T Wnr + r T n T WIIr 
= r T Wr - r T Wnr 


= HW^rlll - 
= l|w 1/2 ©r||2 


r T Wlw T r 

l T w 

|w T r| 2 
l T w ’ 


using the fact that II WFE = WII, thus reproducing ([13- 
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